DOC  FILE  COPY.  ADA038957 


MRC  Technical  Summary  Report  «l7Z-4 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 


610  Walnut  Street 
Madison,  Wisconsin  53706 


ebruary  1977 


^ - 


Approved  for  public  release 
Distribution  unlimited 


and  National  Science  Foundation 

Washington,  D.  C.  20‘iSO 


UNIVERSITY  OF  WISCONSIN  - MADISON 
MATHEMATICS  RESEARCH  CENTER 


ON  MODELUNG  PATTERN  FORMATION 
BY  ACTIVATOR-INHIBITOR  SYSTEMS 

Paul  C.  Fife 


Technical  Summary  Report  «17Z4 
February  1977 

ABSTRACT 

The  formation  of  spatially  patterned  structures  in  biological 
organisms  has  been  modelled  in  recent  years  by  various  mechanisms, 
including  pairs  .of  reaction-diffusion  equations 
u^  = V^u  4-  f(u,v)  , 

Vj.  = V^v  + g(u,v). 

Their  analysis  has  been  by  computer  simulation.  In  some  cases,  u can 
be  interpreted  as  an  activator  and  v an  inhibitor.  The  following  problem 
is  treated:  given  a "pattern"  u = (p{x)  v = i|j(x),  find  a system  which 
has  it  as  a stable  stationary  solution  (stability  is  used  in  various  senses 
in  the  paper).  This  inverse  problem  is  shown  to  have  solutions  for  reason- 


able <p  and  The  solutions  constructed  are  of  activator-in 

with  > D^  . 
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ON  MODELLING  PATTERN  FORMATION 
BY  ACTIVATOR-INHIBITOR  SYSTEMS 

Paul  C,  Fife 

I.  Introduction. 

The  formation  of  spatially  patterned  structures  in  biological  organ- 
isms has  been  modelled  by  various  mechanisms  in  recent  years,  (see, 
for  example,  the  survey  by  Cooke  [2]  for  a description  of  recent  work  in 
the  area).  Among  these  modelling  attempts  have  been  those  using  only  the 
processes  of  chemical  reaction  and  the  diffusion  of  the  reacting  species. 

The  most  common  approach  to  this  problem  has  been  that  of  small-amplitude 
(linear  or  nonlinear)  analysis  of  the  onset  of  symmetry -breaking  instabilities. 

I shall  not  attempt  to  survey  results  using  this  approach,  except  to  say 
that  they  were  begun  by  the  well-known  work  of  Turing,  and  include,  among 
others,  the  work  of  Gmitro,  Othmer,  and  Scriven,  Prigogine,  Lefever,  and 
Nicolis,  and  Segal,  Jackson,  and  Levin.  We  are  concerned  here  with  models 
yielding  large-amplitude  patterns.  Notable  in  this  regard  is  the  work  of 
Gierer  and  Meinhardt  [6,7,9,10]  in  Tubingen  , and  Babloyantz,  Hiernaux, 
Herschkowitz-Kaufman,  Nicolis,  Prigogine,  and  others  in  Brussels  [1,8]. 

See  also  [4]  for  models  of  sharply  differentiated  structures.  These  models 
are  generally  of  the  form 

(la)  u^  = Dj  V'^u  + f(u,v)  , 

(lb)  v^  = V^v  + g(u,v), 
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where  u and  v are  the  concentrations  of  two  hypothetical  reacting 
substances,  and  the  source  terms  f and  g are  derived  from  assumed 
reaction  schemes.  The  analysis  of  these  equations  has  been  by  computer 
simulation  (except  in  [4]).  The  results  have  indeed  indicated  the  appear- 
ance of  spatial  patterns  of  various  sorts. 

The  important  problem  here  is,  given  a stationary  spatial  pattern  of 
some  general  description,  to  find  a system  of  reaction-diffusion  equations 
which  will  yield  that  pattern  as  a stable  solution.  In  addition,  one  wants 
to  be  able  to  interpret  the  reaction  terms  on  the  basis  of  some  reasonable 
reaction  kinetics.  It  is  this  first  inverse  problem  (given  a solution,  find 
the  equations)  that  we  are  concerned  about  in  the  present  paper.  We  re- 
strict attention  to  two-component  systems  of  the  form  (1)  with  one  space 
variable  x (see  §4  for  extensions  to  higher  dimensions),  and  show  that 
for  any  reasonable  given  time -independent  function  u = <p(x),  there  are 
systems  of  this  type  which  have  a stable  stationary  solution  with  u = <p. 
Within  certain  limitations,  the  function  v = 4i(x)  can  also  be  prescribed. 
Systems  having  the  given  functions  as  solution  are  easy  to  construct,  but 
stability  is  an  elusive  property,  and  the  difficult  part  is  finding  systems 
for  which  the  given  solution  is  stable.  For  example,  the  corresponding 
inverse  problem  for  a scalar  equation,  u^  - + f(u),  the  pattern  being 

defined  for  all  x , does  not  have  a solution  unless  the  required 
pattern  is  monotone;  and  even  then  it  is  structurally  unstable, 
as  small  changes  in  f will  destroy  the  pattern  [5]. 
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We  leave  aside  the  second  important  task  mentioned  above,  namely 
the  interpretation  of  f and  g on  the  basis  of  some  reasonable  reaction 
kinetics.  The  systems  we  construct  have  no  obvious  interpretation  in 
those  terms,  and  for  this  reason  are  not  likely  to  be  of  practical  impor- 
tance as  models.  At  the  same  time,  in  constructing  model  reaction-diffu- 
sion systems  with  only  two  components,  one  should  not  attach  overriding 
importance  to  having  them  mirror  specific  reaction  networks  involving  the 
two  species  alone.  In  fact,  the  actual  mechanism  modelled  will  involve 
a large  num.ber  of  reacting  species  (and  other  entirely  different  processes 
as  well).  One  usually  justifies  the  reduction  to  two  species  -ou  the  basis 
that  various  pseudo-steady-state,  slow  reaction,  or  other  approximations 
can  often  be  made  to  effect  such  a reduction.  But  when  these  approxima- 
tions are  introduced,  the  connection  between  the  source  terms  f and  g 
and  the  actual  kinetics  is  necessarily  obscured;  in  particular,  f and  g 
do  not  conform  to  mass -action  kinetics  between  the  two  hypothetical 
species. 

The  only  thing  we  require  of  f and  g is  that  they  be  of  activator - 
inhibitor  type,  which  we  define  as  follows: 

Definition:  u is  an  activator  for  (1)  if  f >0,  g >0. 

V is  an  inhibitor  for  (1)  if  f <0,  g <0. 

Thus  increasing  the  amount  of  u present  enhances  the  production  of  u 
and  V , whereas  increasing  v has  the  opposite  effect.  This  requirement 
was  occasioned  by  the  fact  that  Gierer  and  Meinhardt's  models  are  of 
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activator-inhibitor  type,  with  the  inhibitor  diffusing  more  rapidly  than  the 
activator.  This  latter  fact  is  also  true  in  our  scheme,  and  is  necessary 
to  produce  stability  of  the  pattern.  In  fact,  one  major  point  of  the  present 
paper  is  an  elucidation  of  how  activator-inhibitor  mechanisms  with  differ- 
ent diffusivities  may  enhance  stability. 

Up  to  this  point,  we  have  not  spoken  of  boundary  conditions  which 
one  must  impose  if  (1)  is  to  be  solved  in  a bounded  domain  in  space.  They 
are  of  little  importance  in  our  argument,  which  is  valid  for  any  reasonable 
boundary  conditions,  and  even  for  problems  on  the  whole  line  with  no 
boundary. 

The  concept  of  stability  we  use  for  most  of  the  results  is  linearized 
stability,  wherein  one  examines  the  spectrum  of  the  operator  S obtained 
by  linearizing  the  right  hand  side  of  (1)  about  the  stationary  solution 
If  the  spectrum  S(S)  has  points  \ with  Re  X.  > 0 , then  is  un- 

stable, whereas  if  all  points  have  Re  X.  < 0,  it  is  stable.  Marginal  sta- 
bility is  when  £(S)  has  a point  with  Re  \ = 0,  but  none  with  Re  X.  > 0. 
From  linear  stability,  one  can  often  deduce  stability  in  a stricter  (such  as 
C^)  sense;  in  particular,  that  can  be  done  with  our  examples  of  patterns 
on  a finite  interval,  and  for  our  "single-peak"  pattern  on  the  whole  real 
line.  For  any  pattern  on  the  whole  line,  0 is  always  in  2(S),  so  the 
most  we  can  hope  for  in  the  general  case  is  marginal  stability,  which  is 
what  we  in  fact  obtain.  In  the  case  of  "single  peak"  distributions  on  the 
whole  line,  however,  fo.*  which  <p  is  even  and  monotone  for  x > 0 , it 
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turns  out  that  0 will  be  an  isolated  point  of  ^{S),  and  we  apply  a theorem 
of  Sattinger  [11]  to  obtain  the  result  that  the  stationary  solution  is  stable  in 
the  sense.  I believe  that  our  results  are  among  the  first  in  which 
stability  is  proved  for  solutions  of  systems  with  two  or  more  components, 
on  the  entire  real  line,  or  of  large  amplitude  patterns  on  a finite  interval. 


Some  of  the  results  in  this  paper  were  announced  in  [5]. 

2.  The  inverse  problem  on  a finite  interval. 

We  consider  systems  of  the  form 

(2a)  u^  " ^xx  ^ f(u,v), 

(2b)  Vj.  = k v^  + g(u,v) 

for  functions  u,v  defined  for  xt  [a,b],  t > 0.  We  also  prescribe  bound- 
ary conditions  at  the  endpoints  a and  b , and  for  definiteness  take  them 
to  be  of  no -flux  type: 

(3)  u (a,t)  = u (b,t)  = V (a,t)  = v (b,t)  = 0 . 

X X X X 

For  a given  stationary  solution  u = <^(x),  v = 4/(x)  of  (2),  (3),  we 
denote  by  S the  operator  obtained  by  linearizing  the  right  side  of  (2) 
about  («>,4j). 


^xx  ^ y«i0»4j)v 


kVxx  + g^(<p,4^)u  + g^((p,^)v 


The  operator  S is  to  be  considered  as  acting  on  (C^)^  = C^[a,b]X  C^[a,b], 
2 

with  domain  the  C functions  satisfying  ( 3).  We  denote  its  spectrum  by  Z(S). 

2 

The  patterns  <p{s)  which  we  consider  will  be  in  C [a,b],  and 
will  have  the  property  that  </>"  is  a function  only  of  tp  . This 
latter  means  that  ^ is  even  with  respect  to  any  local  maximum  or 
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minimum,  so  that  its  graph  (see  Figure  1)  is  symmetric  with  respect  to 
vertical  lines  through  such  maxima  and  minima.  To  be  more  precise,  the 
statement  is  that  <p  can  be  extended  to  be  defined  for  all  x , and  the 
extended  graph  has  those  symmetry  properties. 


Figure  1.  Example  of  a pattern  ^(x)  with  two  vertical  axes  of 

symmetry  (dotted  lines). 

' Z 

Theorem  1;  Let  ^ be  a function  in  C [a,b],  satisfying  <p'(a)  = <p'{h)  = 0, 

2 

and  <p"(x)  = F{(p  (x))  for  some  C function  F . Let  satisfy 
(4)  <p  = Cj4)  + c^  , c^  > 0 . 

Let  k > 1 . 

Then  there  exist  functions  f(u,v),  g(u,v),  satisfying  f^  >0, 

>0,  < 0,  g^  < 0,  such  that  is  a stable  or  marginally  stable 

stationary  solution  of  (2),  (3).  Furthermore  2(S)  is  discrete  and  is  located 
entirely  in  the  half-plane  {Re  \ < 0}  except  possibly  for  a simple  eigenvalue 
at  the  origin. 
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Proof : First  consider  the  case  so  c^  = 1,  c^  = 0 . For  some 

constant  c to  be  determined  later,  let 
(5a)  f(u,v)  = -F(u)  + <r(u-v) 

(5b)  g(u,v)  = -k  F(v)  + ko-(u-v)  . 

If  (T  is  chosen  large  enough,  these  will  be  the  required  functions.  First  of 
all,  it  is  clear  that  (u,v)  = (<f>,<p)  satisfies  (2),  (3).  We  must  show  it  is 

d2 

a stable  or  marginally  stable  solution.  Let  L = — •_  - F'(v'(x)),  an 

0 2 
operator  on  C [a,b]  with  domain  fi(L)  consisting  of  C functions 

satisfying  (3).  Then 


/ L u + o-(u-v)  \ 


ykL  V + ko-(u-v) y 
Let  e = l/k  , and  let  P(X.,(j.)  be  the  polynomial 

P(k,p.)  = e + \[(l-e)o-  - (l+e)p.]  + . 

Then  P(X,  L)  is  a fourth  order  operator  with  domain  functions  he  fi(L) 

such  that  Lh  e fi(L);  furthermore  S(P(X.,L))  = P(\,Z(L))  [3,  p.  604].  Let 

A = Z:(P(\,L))  9 0}  = {X  : P(X,p)  = 0 for  some  p e 2(L)  }. 

Proposition  1 ; Z(S)  C A . 

Proof:  Let  X ^ A . Then  0 is  in  the  resolvent  set  of  P(X,L),  and 

P(X,L)  ^ is  a bounded  operator  on  C^[a,b].  Consider  the  problem 


(6) 


= 1^1, 


with  p,q  e C^[a,b],  It  can  be  solved  explicitly  by  the  formula 
(7a) 

(7  b) 


u = (L  - <r  - eX)P(X,L)'^p  + e 0-  P(X,L)‘^q  , 


V = -0-  P(X,L)‘V  + e(L  + o--X)P().,L)‘^q  , 
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commutes  with  L on 


It  is  verified  directly,  using  the  fact  that  P(\,L) 
fi(L),  that  this  formula  defines  a two-sided  bounded  inverse  (S-X)  There- 
fore X is  in  the  resolvent  set  of  S , This  establishes  the  proposition. 
Proposition  2;  A is  a discrete  set  in  the  half-plane  (ReX  <0}.  Further- 
more A C {ReX  <0}  unless  0 « 2(L),  in  which  case  A = {0}  U T , 
with  T C (ReX  < 0 }. 

Proof:  The  discreteness  follows  from  the  discreteness  of  2(L).  We  know 

that  21(L)  is  real  and  bounded  from  above.  Let  a be  so  large  that 

(l-e)(r  - (l+e)Z(L)  > 0.  Then  if  p.  t 2(L)  and  p 0,  the  coefficient  of  X 

in  the  polynomial  P(X,p)  will  be  positive,  so  that  the  sum  of  the  two  roots 

^I’^Z  ~ ® negative.  If  they  are  complex  conjugates, 

their  common  real  part  must  be  negative.  If  they  are  real,  then  X^X^  = 

2 

p > 0,  so  that  they  must  both  be  negative.  In  any  case.  Re  X < 0 for  any 
root  of  P . 

On  the  other  hand  if  p = 0 e 2(L),  then  P(X,0)  has  one  negative 
and  one  0 root.  This  proves  the  proposition. 

We  now  observe  that  f^  > 0,  etc.  , for  large  enough  o-  . To  com- 
plete the  proof  of  Theorem  1,  we  need  only  show  that  when  0 is  an  eigen- 
value of  S , it  ^s  simple.  From  Proposition  2,  we  know  this  happens  only 
when  0 is  an  eigenvalue  of  L . Under  boundary  conditions  (3),  L has 
only  algebraically  simple  eigenvalues  so  there  is  a unique  (up  to  a scalar 
factor)  nullvector  0 . Let  (Xj,X2)  be  any  nullvector  of  S , so  that 
LXj  + o-(Xj-X2)  = 0 > 

9 

LX2  + <r(X2-X2)  = 0. 
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Subtractinq,  we  see  that  L(\  -x  ) = 0 , which  implies  that  x -X  = a0  for 
some  scalar  a.  Hence  Lx^  + aaO  = 0.  But  since  0 is  an  algebraically 
simple  eigenvalue  of  L , this  implies  a = 0 and  = bO.  Hence 
(XjtX^)  = b(G,0),  and  S has  only  one  nullvector.  A similar  argument  shows 
that  its  algebraic  multiplicity  is  also  1 . This  completes  the  proof  of  the 
theorem  in  the  case  (;(’  = 4j. 

Now  if  <p  * i\>,  but  is  given  by  (4),  we  first  write  the  system  (2) 
with  f and  g given  by  (5),  and  the  symbol  v replaced  by  w.  Then  we 
effect  a change  of  variable  w = c^v  + c^  to  obtain  the  desired  system  in 
u and  V . This  completes  the  proof. 

Remarks : 1.  If,  instead  of  (4),  tp  and  4j  are  related  nonlinearly  by 

<p  = h(vl;) 

with  h'(4j)  > 0,  then  we  still  obtain  a type  of  nonlinear  diffusion  system  for 
u and  V.  The  difference  is  that  now  (2b)  is  replaced  by 

(h'(v)Vx)x 

''t ' — i?(7j — ♦ • 

This  is  obtained  by  the  procedure  discussed  in  the  preceding  paragraph  , 
except  that  the  change  of  variable  is  now  w = h(v). 

2.  The  specific  form  of  the  boundary  conditions  (3)  was  not  used; 
only  the  fact  that  21(L)  is  bounded  from  above,  discrete,  and  contains  0 , 
if  at  all,  only  as  a simple  eigenvalue.  If  this  last  condition  is  not  fulfilled, 
then  all  assertions  of  the  theorem  remain  true  except  the  final  one  regarding 
the  simplicity  of  the  eigenvalue  at  the  origin.  We  may  therefore  replace  (3) 
by  any  other  bourfdary  conditions  for  which  the  above  requirements  are  met. 
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In  particular,  periodic  boundary  conditions  could  be  used.  In  this  case,  0 
is  always  an  eigenvalue  of  L , because  L(p'  =0.  If  L has  constant 
coefficients,  it  would  of  course  be  double,  because  L<p"  = 0 as  well.  But 
in  general,  one  expects  it  to  be  simple. 

3.  The  solution  of  the  inverse  problem  is  certainly  not  unique;  for 
example,  o-  can  be  any  large  enough  positive  constant.  More  generally, 
in  place  of  (5),  we  may  use  functions 
(8a)  f(u,v)  = -F(u)  + (t(u-v)  + (u-v)^  f(u,v) 

(8b)  g(u,v)  = -lcF(v)  + Ict(u-v)  + (u-v)^  g(u,v), 

where  f and  g are  arbitrary,  o-  and  t are  such  that  there  ^re  no  points 
of  2(L)  in  [T-o-,0)(if  (t>t)  or  in  (0,t-(t]  (if  er  < t),  and  t > e cr  + (1+e )M, 
where  M = Max(2(L)).  With  such  a pair  f, g,  the  proof  of  Theorem  1 goes 
through  with  only  obvious  modifications. 

3.  Patterns  on  the  entire  line. 

The  restriction  to  a finite  interval  in  the  previous  section  is  un- 
necessary. Let  (|p  be  a bounded  nonconstant  function,  defined  for  all  x , 

2 

such  that  ip"  = F(<p)  for  some  C function  F . Then  (p  will  be  periodic, 
pealced  or  monotone.  Here  "peaked"  means  that  it  has  a single  maximum  or 
minimum,  and  approaches  a limit  as  |x|-»  oo  . 

Theorem  2.  Let  <p  be  as  described,  and  let  ijj  satisfy 

<p  = c^ijj  + c^,  Cj^  > 0 . 

Let  k > 1,  Then  the  conclusions  of  Theorem  1 hold,  except  that  2(S)  is 
no  longer  discrete,  and  in  the  periodic  case,  the  eigenvalue  at  the  origin 
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is  not  necessarily  simple.  In  all  cases,  0 c 2(S).  In  the  peaked  case,  let 
<p  = lim  <p(x).  If  F'(</J.)  > 0 , then  0 is  an  isolated  point  of  i;(S)  and 
{<p,  4j)  is  steLle  in  the  uniform  norm.  In  the  monotone  case,  let  = 
lim  v>(x)  and  = lim<;()(x).  If  F'((j()^)  > 0 , i =0,1,  then  is 

X-*  -00  X-*-  00 

again  stable  in  the  uniform  norm. 

Proof ; The  proof  of  Theorem  1 holds  without  change,  with  the  same 
functions  (5),  except  that  A will  no  longer  be  discrete,  because  2:(L) 
is  no  longer  discrete.  In  the  peaked  case  with  F'(<^q)  > 0 , we  know  that 
the  points  of  Z(L)  greater  than  -F'(<;«>q)  are  isolated.  This  follows  from 
[3,  p.  1448]  and  the  easily  proved  fact  that  any  bounded  solution  of 
(L-\)u  = 0 for  X > -F'(v’q)  must  decay  exponentially,  so  is  in  y^(R): 
hence  discrete  points  of  L in  this  range  are  eigenvalues  of  L as  an 
operate  in  other  hand,  0 is  always  an 

eigenvalue  of  L,  since  L<p'  = 0 , and  it  is  simple  and  has  a finite  number 
of  nodes.  So  although  A is  no  longer  discrete,  it  still  lies  in  the  left  half- 
plane except  for  an  isolated  point  at  the  origin.  The  same  is  therefore  true 
of  ^{S).  The  proof  that  0 is  a simple  eigenvalue  of  2(S)  proceeds  as 
before. 

To  prove  that  is  stable  in  the  sense,  we  shall  use 

Sattinger's  stability  theorem  [11].  To  verify  that  theorem's  hypotheses, 
it  suffices  to  show  that  i:(S)  lies  on  the  negative  real  axis,  and  that  the 
resolvent  (S-X)"^  satisfies  the  following  estimates  for  ] arg  Xj  <tt-6  < it 
and  |x|  large  enough.  (In  the  following,  the  symbol  K will  denote 
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several  different  constants. ) 


(9a) 

(9b) 


ll(S-M'^h||  < 

(C  ) 

ll(S-M’^h||  < 

(C')^ 


M IlHI 


0 2 

(C  ) 


0 2 

(C  ) 


> 


It  may  be  verified  directly  from  the  definition  of  A,  and  the  fact 
that  2(L)  is  real,  that  A,  hence  2(S),  is  real  for  o-(l-e)  large  enough. 
Therefore  2:(S)  certainly  lies  in  the  required  sector. 

We  now  turn  to  the  verification  of  (9).  First,  we  note  that  esti- 
mates of  that  form  were  obtained  in  [12]  for  the  operator  L instead  of  S , 
provided  F'(^(x))  approaches  its  limit  exponentially  fast  as  |x|-*oo.  But 
F(<Pq)  = 0 ) ^'(<^0^  ^ satisfies  <p"  - F(<^)  = 0 , so  «p(x)-» 

exponentially  as  |x|-*oo.  Therefore  F'(<p(x)) -►  F'(<Pq)  exponentially  as 
well . 


1 

1?1 


We  may  factor  P(\,L)  as  follows: 

P(\,L)  = (L-Hi^KL-Hi^)’  = i[Ml+e)±  [ V^(l-e)^-4<r(l-e)\]']. 
It  is  seen  that  for  large  | X| , are  also  large  and  for  some  K , 


(10) 


.-1 


K'  IM  < IpJ  < K|xl  . 


-1 


By  (9)  and  the  above  mentioned  estimate  for  (L-\)  , we  have  the  following, 

for  large  enough  ( X ( in  the  sector  ( arg  X | < ir  - 6 : 

K 


(11) 


If  (L-^^)w  = h,  we  have  w"  = h - p.^w  + F'(^(x))w  , so  from  (11)  , 

2 - 0 

’ ^ C^(R)  C (R) 
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Now 


lP(\,L)'^h| 


c 


and  so 


< K||(W^)-'h||  „ < i^llhll  „ , 


|LP(\,L)'^h|l  Q < K|lP(\,L)‘^h|l  2 < 


These  estimates,  applied  to  the  explicit  representation  (7)  of  (S-X)  \ 
yield 

(U)  ■ 

up  - 

If  (S-X)(  - ) = ( ),  we  may  again  represent  u”  and  v in  terms  of 

V q 

u,  V,  p,  and  q,  and  so  obtain 

Interpolating  between  this  estimate  and  (12),  we  finally  obtain 

This  establishes  the  C^-stability  of  peaked  distributions. 

For  monotone  distributions,  S(L)  is  discrete  above  Max[-F'(v’Q)^ 
-F'(^^)]  < 0 , The  rest  of  the  argument  is  the  same.  This  completes  the 
proof  of  Theorem  2. 

Example ; Possibly  the  simplest  illustration  of  single-peak  modelling 
would  be  for  the  pattern 


<p{x)  = 


(1+e  ' ) 


6 e ' 


which  represents  p peak  with  maximum  of  3/2  at  x = 0,  <p{<x>)  = 0 , and 
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- 2 2 

an  approximate  width  d , which  is  arbitrary.  Since  <p"  + d {(p  -<p)  - 0 , 

the  model  reaction-diffusion  equations  can  be  taken  to  be 

-2  2 -2 

U,  = U + d U + ((T-d  )u  - (TV  , 

t XX 

-2  2 -2 

V.  = kv  + kd  V - (ff+kd  )v  + o-u  . 

t XX 

Remark : In  the  peaked  case,  there  is  another  stationary  solution  of  (2) 

with  (5)  besides  u = v (x),  namely  the  constant  solution  u = v = , 

It  is  also  stable,  as  the  above  analysis,  together  with  the  fact  that 
F'(<Pq)>  0 , shows.  This  suggests  there  may  be  a third  stationary  solution 
which  is  unstable.  In  fact,  the  function  F(^p)  must  take  on  positive  and 
negative  values  for  so  there  must  be  another  zero  of  F* , say  <p^  , 

for  which  F'(<p^)  < 0 . Then  u = v = ^ is  a solution.  If  S is  the  lineariza- 
tion about  this  solution,  the  above  analysis  shows  its  spectrum  to  be  con- 
tinuous and  to  extend  up  to  the  origin.  Very  possibly  this  third  solution, 
though  marginally  stable  in  the  linearized  sense,  will  turn  out  to  be  unstable 
in  the  sense. 

4.  Discussion  ; 

I.  Many  of  the  results  extend  to  analogous  problems  in  more  than  one 
dimension.  For  this  extension,  one  replaces  u and  v in  (2)  by  V^u 

’ XX  XX 

2 2 
and  V V respectively,  and  considers  patterns  <p  for  which  V <p  - F(^), 

Single  radially  symmetric  peaks,  for  example,  have  this  property,  and  it  is 

likely  that  lattices  of  peaks  can  also  be  constructed  with  it.  Then  i^(S) 

again  lies  on  the  negative  real  axis.  However,  the  eigenvalue  0 is  no 

longer  necessarily  simple. 
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Z.  It  is  clear  that  the  stability  proofs  given  in  the  preceding  two 
sections  depend  on 

(i)  large  enough  activation  by  u and  inhibition  by  v , and 
(ii)  the  diffusivity  (k)  of  v being  larger  than  that  (1)  of  u . 

At  the  same  time,  it  is  clear  from  Remark  3 at  the  end  of  §2  that  the  param- 
eters O'  and  T cannot  be  completely  independent  for  the  stability  proof 
to  go  through.  We  interpret  this  by  saying  that  the  u-activation  and  the 
v-activation  cannot  be  completely  independent. 

In  our  constructed  models,  there  is  also  a relation  between  the 
magnitudes  of  the  activation  and  the  inhibition,  and  this  relation  depends 
on  the  relative  amplitudes  of  the  required  distributions  q>  and  ijj.  For 
simplicity,  let  us  take  c^  = 0 in  (4);  then  this  relative  amplitude  will  be 
c^.  For  Cj  not  necessarily  1 , the  equations  corresponding  to  (8)  will  be 

f(U,V)  = -F(U)  + (t(U-CjV)  + (u-c^v)^  f 

k kx  Z 

g(u,v)  = --F(c  v)+  — (u-c.v)  + (u-c.v)  g . 
c^  1 c^  1 1 

It  may  be  instructive  to  express  all  the  relations  mentioned  above 
directly  in  terms  of  some  reasonable  activation  and  inhibition  parameters. 
Accordingly,  we  define  the 

u-activation  - = average  value  of  f^(v>,4j)  ~ o-  (for  o-  large), 

kx 

v-activation  = average  value  of  g^(  — (for  r large), 

u-inhibition  = I ~ c.cr 
u 1 

v-inhibition  = I ~ kx  . 

V 
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For  our  models,  then,  the  following  relations  exist  between  the 


parameters  A , A , I , I , c,,  and  k : 

U V u V 1 

A ,A  ,I  ,r  » 1 , 

u v'  U V 

k > 1, 


(13) 


c.A  , 
1 u 


I ~ c.A 
V 1 V 


A ~ — A 


V ^ 

All  of  these  are  clear  from  the  preceding  except  possibly  the  last.  It  follows 
from  the  fact  that  although  a and  t were  required  to  be  large,  their  differ- 
ence had  to  be  small  enough.  So  we  can  write  <r  ~ t . 

This  is  no  demonstration  that  the  relations  (13)  are  always  necessary 
to  produce  stable  patterns;  nevertheless  it  seems  reasonable  to  use  them  as 
one  possible  guide  in  constructing  other  models,  not  of  the  form  given  in 
the  present  paper. 
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